Tunable interactions and phase transitions in Dirac materials in a magnetic field 
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A partially filled Landau level (LL) hosts a variety of correlated states of matter with unique 
properties. The ability to control these phases requires tuning the effective electron interactions 
within a LL, which has been difficult to achieve in GaAs-based structures. Here we consider a 
class of Dirac materials in which the chiral band structure, along with the mass term, give rise to 
a wide tunability of the effective interactions by the magnetic field. This tunability is such that 
different phases can occur in a single LL, and phase transitions between them can be driven in situ. 
The incompressible, Abelian and non-Abelian, liquids are stabilized in new interaction regimes. 
Our study points to a realistic method of controlling the correlated phases and studying the phase 
transitions between them in materials such as graphene, bilayer graphene, and topological insulators. 
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In a magnetic field, the kinetic energy of electrons 
confined to move in two dimensions is quenched into a 
set of discrete Landau levels (LLs). The properties of 
a partially filled LL are therefore determined solely by 
the electron-electron interactions, which give rise to a 
number of fundamentally new many-body phases, includ- 
ing the incompressible fractional quantum Hall (FQH) 
states^— , compressible Fermi-liquid-like states (CFL)£, 
as well as states with broken translational symmetry, 
such as charge-density waves (CDW), stripes and bub- 
ble phases^. 

The competition between different phases at a given 
partial filling v is sensitive to the form of the effective 
Coulomb interactions within a LL. Experimentally, to 
date there is no reliable way of controlling the effective 
interactions within a LLi. This is because in GaAs, the 
most common high-mobility 2D electron system, the ef- 
fective interactions depend only on the LL number, and 
not on the magnetic field. Thus, in order to control the 
FQH phases, it is advantageous to find alternative 2D 
electron systems with tunable interactions^. 

Recently, a new class of such materials clean enough 
to observe FQH phases has emerged. These so-called 
Dirac materials host chiral excitations with the non- 
trivial Berry phases. Examples include graphene, bilayer 
graphene 9 , topological insulators 10 , as well as certain 
quantum wells^i. A natural question arises - could the 
Dirac materials offer any tunability in the FQH regime? 

Here, we answer this question affirmatively, and pro- 
vide a realistic model where the effective interactions can 
be widely tuned by varying an external magnetic field, 
giving rise to several phases within each LL, and quan- 
tum phase transitions (QPTs) between them. We map 
out the phase diagram of chiral materials at certain fill- 
ing factors, identifying a new regime of the effective in- 
teractions where the non- Abelian Moore- Read (MR) 12 
and other paired states^ are stable away from n = 1 
LL of GaAs^r— . Finally, we predict several types of 
phase transitions that occur in Dirac materials (from an 



incompressible FQH state to a compressible state, with 
or without breaking of translational symmetry; between 
an Abelian and a non-Abelian FQH state), which previ- 
ously have only been considered in artificial theoretical 
models. 

The control of the correlated phases proposed here is 
attractive for two reasons. First, it allows one to real- 
ize and stabilize the exotic states^, or may lead to a 
discovery of new correlated states that are weak or ab- 
sent in GaAs. Second, it provides a setting for studying 
the fundamental problem of phase transitions that in- 
volve topologically ordered states. In our model, both 
of these goals can be achieved by current experimental 
techniques of tuning the gap in the spectrum, e.g. in bi- 
layer graphene 9 -, and topological insulators 10 . In the for- 
mer case, the gap is opened by a perpendicular electric 
field, and in the latter case by the deposition of magnetic 
adatoms. In graphene, it is more challenging to open 
the gap, although there exist several promising propos- 
als (e.g., the mass can be generated either spontaneously, 
or as a result of sublattice symmetry breaking 9 -). 

We consider fermions with the Berry phase 7r and 27r, 
described by the 2x2 Hamiltonian, 



A = 1,2 



A Mx(p x +ip y y 
Mx(p x -ip y ) x -A 

(1) 

where A^i = v is the Fermi velocity, A4 2 = l/2m (m 
is the effective mass), and 2 A is the band gap. The case 
of 7r-carriers (A = 1) is realized in graphene, topological 
insulators, and special quantum wells-ii; the case of 27r- 
carriers (A = 2) occurs in bilayer graphene^. 

Landau level spectrum for 7r-carriers, obtained by solv- 
ing the Schrodinger equation for H n in a magnetic field, 
is given by e„ = sgn(n)^/A 2 + e§|n|, n = ±1,±2,..., 
£o = V2hvo/ (-B is the characteristic energy scale, and 
£b = y/Hc/eB is the magnetic length. The correspond- 
ing two-component wave functions are given by ip n — 
(cos9 n ipi n \_x, shid n <p\ n \), where ip n is the wave function 



2 



1 

0.8 
a. 0.6 



CD 



0.4 



0.2 



°N=3 ■ Y 




n N=9 




-°N=10 - 




A N=8 - 




A N=9 | 




- A N=10 "" 




- (st,2) 

f | : 'l 

I P ( 1 


(Jt,l) 


[■'!■ J 





0.1 

0.08 

0.06 



0.04 ffl 



0.0 0.5 1.0|0.0 0.5 1.0 



0.02 




FIG. 1. (Color online) Transition between the Laughlin and 
the bubble state. Overlap with the Laughlin wavefunction 
On and charge gap Ajv are plotted for an iV-particle system 
at v = 1/3 filling of (n, 1) and (tt, 2) LLs. 
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FIG. 2. (Color online) Transition between the Moore-Read 
and the CFL state. Overlap with the Pfaffian wavefunction 
On and charge gap A at are plotted for an iV-particle system 
at v = 1/2 filling of (tt, 1) LL. 



of the nth non-relativistic LL (standard magnetic oscilla- 
tor wave function) , and parameter 9 can be expressed as 

tan0 n = sgn(n) v /(A/e ) 2 + \n\ - A/e Q / y/\n\. In the 
limit of zero mass (graphene case), tant9„ = ±1, and the 
weights of V|n|-i an d Vlnl m the wave function ip n become 
equal. In the opposite limit of very large mass, A/eo 3> 1, 
the LLs become identical to the non-relativistic ones: 
tan#i n i — > (this corresponds to tp\ n \ « (y>i n i_i , 0)), and 
tant9_i n i — > — oo (V*— inl ~ (0) f\n\))- As we comment 
below, varying A/eo between and oo, which can be 
achieved by changing magnetic field, allows one to ex- 
plore the whole range 9 6 (0; 7r/2). 

The two-component nature of the wavefunction mod- 
ifies the effective interaction within a LL, which deter- 
mines the many-body phases at a partial filling. We 
use the standard approximation and project the inter- 
action onto a partially filled LL. In this case, the Fourier 
transform of the effective interaction is a product of bare 
Coulomb interaction, V(q) = 27re 2 /g, and the form fac- 
tor— |F n (<y)| 2 which contains the information about the 
band structure. In case of the materials with the Berry 
phase tt, 

F£(q) = cos 2 9L lnhl (q 2 l%/2) + sin 2 9L H (q 2 £ 2 B /2), (2) 

where Lk is the fcth Laguerre polynomial, and for sim- 
plicity we omitted the index of 9. The form-factor is 
a mixture of the (|n| — l)th and |n|th LL form-factors 
in a non-relativistic 2DES with parabolic dispersion. At 
9 = 7r/4, the above equation reduces to the form- factor 
of pristine graphene^. 

Similarly, for carriers with the Berry phase 27r, the LL 
spectrum is given by e„ = sgn(n) -\/ A 2 + e;?|ra|(|n| + 1), 
n = ±1,±2,..., where s c — eB/mc is the cyclotron 
energy. The corresponding wave functions are 
ip n = (cos6 n ip\ n \_i,sin9 n ip\ n \ + i), with tan#„ = 

"sgn(n) V(A/e c ) 2 + NCR + 1) " A/e c ] / y/\n\(\n\ + 1). 
The form-factor is a mixture of standard (|n| — l)th and 



(|n| + l)th form-factors, 



F 2 *(q) 



cos 



] 6L lnhl (q 2 e%/2) 



i 2 9L H+1 (q 2 e 2 



(3) 

The tunable form of the effective interactions Eqs. (l2l3[ ) 
provides a way to analyze the transitions between 
strongly correlated phases in a direct manner, as will be 
demonstrated below. 

Much of the previous theoretical work on the quantum 
Hall effect in the Dirac materials has been limited to 
graphene, exploring the consequences of the four-fold LL 
degeneracy (valley and spin) that leads to new SU(2) 
and SU(4)-symmetric fractional states^. Here we neglect 
the multicomponent degrees of freedomi^, and examine 
the effects originating from the interplay of the Coulomb 
interaction and band structure. The large variation of 
the effective interactions, due to the band structure, is 
assumed to be the dominant effect, even when corrections 
due to LL mixingi^ are taken into account. 

We proceed by studying the interacting states of Dirac 
materials using exact diagonalization in the spherical and 
torus geometr y 20 i 21 . The former is useful in studying the 
incompressible liquids, but is not suitable for states that 
break translational symmetry, when periodic boundary 
conditions are more natural 2 ^. If the system has an un- 
derlying lattice structure, it is assumed that its lattice 
constant is much smaller than is- We map out the phase 
diagram as a function of 9 using the overlaps between 
an exact ground state and a trial wavefunction. More- 
over, we evaluate the energy gaps for creating charged 
excitations, A = E qp + E q h — 2E , where E q h(E qp ) is 
the energy of the system in the presence of a quasihole 
(quasiparticle) and Eq is the ground state energy (all in 
units of e 2 /e£b)- For simplicity, when we use the spher- 
ical geometry, the form factors Eqs.(0 [3]) are calculated 
for a flat surface; the curvature is expected to produce 
only small quantitative difference in larger systems and 
has a stronger impact on the energy gaps than the over- 
laps. Various choices of the single-particle Hamiltonian 
are denoted by (Xtt, n), where we restrict to A = 1,2 and 
n = 1,2. Higher values of n are not considered because 
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FIG. 3. (Color online) Phase diagram of a generic chiral material in a partially filled LL containing a mixture of n = 0, n = 1 
and n = 2 non-relativistic form-factors. Color scheme on the right defines the magnitude of the overlaps with the Laughlin, 
Moore-Read and the k = 3 Read-Rezayi (RR) wavefunction. Shaded regions represent the phases with different topological 
numbers (CFL, hierarchy state). The phase boundary between topological states is approximately drawn by a wiggly line. 



they do not support FQH states. 

At partial filling v — 1/3, the Laughlin state is robust 
in (tt, 1), (tt, 2), and (2tt, 1) LLs. The overlaps between 
the exact ground state and the Laughlin wavefunction, 
as well as the charge gaps in (71", 1) and (tt, 2) LLs, are 
illustrated in Fig. [TJ The overlap remains very close to 
one for cos 2 8 £ [0.2; 1] in the (tt, 1) LL. This is not sur- 
prising because the form-factors of this LL are a mixture 
of and 1 non-relativistic form-factors, both of which fa- 
vor the Laughlin state. An interesting feature of Fig. [T]is 
that the energy gap shows a maximum at cos 2 8 rs 0.65, 
rather than at cos 2 8 = 1 (pure n = non-relativistic 
LL). This special value of 8 maximizes the ratio of V1/V3 
Haldane pseudopotentials- which controls the gap in this 
case. The overlap and the gap decrease drastically as 
cos 2 8 — > 0; furthermore, in (tt, 2) LL, the Laughlin state 
quickly undergoes a QPT to a bubble phase as the over- 
laps and gaps drop to zero. A signature of this transition 
is also detected in the projected structure factor—, which 
develops a sharp peak in the bubble phase 2 ^. The tran- 
sition is more naturally captured in the torus geometry, 
where a large ground state degeneracy 2 ^, characterized 
by a 2D array of crystal momenta, develops in (tt, 2) 
LL. Similar behavior is found in (27r, 1) LL 2 — . One of 
the of the members of the multiplet belongs to the zero- 
momentum sector, which suggests that the transition is 
likely second-order. 

Next, we consider a half-filled LL, where we find evi- 
dence for the Moore-Read^ correlations in half-filled chi- 
ral (tt, 1) and (27T, 1) LLs. In the former case, the Pfaf- 
fian becomes more robust than its GaAs analogue for 
cos 2 8 ~ 0.1, Fig. [3J The point cos 8 = corresponds 
to a n = 1 non-relativistic LL where the overlap is high 
(as expected). As cos 6* increases, the overlap reaches a 
maximum value at cos 2 8 m 0.1; this is followed by an 
enhancement of the gap, which peaks at approximately 
the same point (gaps show stronger finite-size effects than 
the overlaps). The increased stability of the Pfaffian was 
previously discussed in Ref.— in the context of biased 
bilayer graphene. On the other hand, away from the 
"optimal" point cos 2 8 0.1, the Pfaffian undergoes a 
QPT to the compressible CFL states*^. This transition 



is accompanied by a change in shift 25 of the ground state 
beyond cos 2 8 = 0.4, and is also manifested in the gaps 
dropping to zero (FigJ^J). A similar scenario is found in 
(2-7T, 1) LL3. At larger values of cos 2 6*, in both (tt, 1) 
and (2-7T, 1), we expect a QPT into the CFL state; at 
very small cos 2 8 in (2-7T, 1) LL, the Pfaffian gives way to 
a stripe phase. Thus, chiral materials are suitable for 
studying phase transitions between the Pfaffian, stripe 
phase and the CFL; these transitions are expected to be 
of different nature, either first or second order—. This is 
a distinct advantage over GaAs-based 2DES where the ef- 
fective interactions are significantly more difficult to tune 
in a controlled fashion. 

Note that the enhancement of the overlaps and gaps 
might appear similar in nature to tweaking of the V\ 
pseudopotential, known to have a favorable effect on the 
Pfaffian stated. However, in the present case, the in- 
teraction that favors the MR state is a superposition 
of n = and n = 1 form factors, including the long 
range tail of the repulsive potential; thus it represents 
a new regime where the MR state is stable. This new 
regime of stability does not crucially depend on the pres- 
ence of n = 1 LL form-factor. Similar signatures of the 
paired states occur also in (27T, 1) LL, where the inter- 
action involves a superposition of n = and n — 2 LL 
form-factors, making the difference from GaAs even more 
striking. 

We have also analyzed the filling v = 3/5^2 (or v = 
2/5), where the one expects to find the non-Abelian k = 3 
Read-Rezayi stated that supports universal topological 
quantum computation 2 ^. Similar to the MR state, in 
chiral materials we find non-Abelian correlations in both 
(tt, 1) and (2-7T, 1) LLs, with a phase transition to the 
stripe phase and the Abelian hierarchy stated (see also 
Fig®. 

Finally, we consider a model that generalizes the effec- 
tive interaction given above and allows to map the phase 
diagram for a wider class of materials. Specifically, we 
study the form-factors that are represented as a linear 
combination of \n\, \n\ + 1 and |n| + 2 non-relativistic 
form-factors. Such form-factors define a two-parameter 
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family, 

F n (q) = cos 2 6L H (q 2 £ 2 B /2) + sin 2 9 cos 2 <pL H+1 (q 2 £ 2 B /2) 
+ sin 2 0sm 2 ^|„| +2 (g 2 4/2). (4) 

The effective interactions of the above form arise in a 
number of materials, including trilayer graphene 2 ^, as 
well as bilayer graphene in the limit of large asymme- 
try between the two layer o 24 ' 29 . The exact relation of 
the parameters (9, ip) to the band structure of various 
materials will be discussed elsewhere^ 3 -. The phase dia- 
gram of the model in Eq. Q for its n = LL is illustrated 
in Fig. [3] Along certain axes (indicated by arrows) , the 
generalized model reduces to one of the particular cases 
(Att, n) presented earlier. A salient feature of the general- 
ized model is that non-Abelian states are found in a strip 
of 9, ip values where the effective interaction significantly 
deviates from n = 1 non-relativistic one. On the right 
side of the strip, the states undergo a QPT to another 
FQH (or CFL) state with different topological numbers. 
These transitions are represented by wiggly lines in Fig. 
121 where the shaded areas mark the resulting phases with 
a different shift. On the left side, the system crosses over 
to compressible, CDW-like phases. The Abelian states 
dominate over a wide region of parameter space and are 
insensitive to variation in ip, unless cos 9 is close to zero. 
We note that the presented phase diagram is consistent 
with the one obtained by considering the charge gaps 



instead of the overlaps 23 . 

In summary, we showed that the new regimes of the 
effective Coulomb interactions can be realized in chiral 
materials. This allows one to stabilize the desired phases 
(including non-Abelian ones) within a single LL, and pro- 
vides a way to engineer QPTs between them. Our results 
apply to a number of available high-mobility 2DES^ i 10 ' 28 , 
however bilayer and trilayer graphene appear as most 
suitable candidates because their band structure can be 
tuned by an external electric field. In fact, our pro- 
posal can be already be realized in bilayer graphene with 
the existing experimental tools. It was demonstrated 3 - 3 * 
that the gap can be controlled in the interval ±125meV. 
With the cyclotron energy of ~ 20meV at B = 10T 
this translates into the ability to tune cos 2 9 in the inter- 
val (0.01; 0.99), sufficient for the realization of the non- 
Abelian states (see Fig. [3]). We believe that experimen- 
tal observation of the phases and QPTs predicted here 
should be feasible in the near future, as the fractional 
quantum Hall states in graphene have already been ob- 
served 3 ^, and transport anomalies suggesting fractional 
states have been seen in topological insulators^! and bi- 
layer graphene 3 - 2 -. 
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